arXiv: 1505.04314vl [nlin.CD] 16 May 2015 


arXiv.org (J. Exp. Theor. Phys., 2015, Vol. 120, No. 5, pp. 912-921) 


MSG: 37M10, 37M20, 34C28, 34C15, 34D06; 

PACS: 02.70.-c, 05.45.-a, 05.45.Tp, 05.45.Xt. 

Analysis of the Time Structure of Synchronization 
in Multidimensional Chaotic Systems 

A. V. Makarenko®’ 0 

® Constructive Cybernetics Research Group 
P.O.Box 560, Moscow, 101000 Russia 

^ Institute of Control Sciences, Russian Academy of Sciences 
ul. Profsoyuznaya 65, Moscow, 117977 Russia 

Received October 28, 2014 

Abstract. A new approach is proposed to the integrated analysis of the time struc¬ 
ture of synchronization of multidimensional chaotic systems. The method allows one 
to diagnose and quantitatively evaluate the intermittency characteristics during syn¬ 
chronization of chaotic oscillations in the T-synchronization mode. A system of two 
identical logistic mappings with unidirectional coupling that operate in the developed 
chaos regime is analyzed. It is shown that the widely used approach, in which only 
synchronization patterns are subjected to analysis while desynchronization areas are 
considered as a background signal and removed from analysis, should be considered as 
methodologically incomplete. 
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1. Introduction 

Synchronization is one of the fundamental concepts of the theory of nonlinear dynamics 
and chaos theory. This phenomenon is widespread in nature, science, engineering, and society [1]. 
One of important manifestations of this phenomenon is the synchronization of chaotic oscillations, 
which was experimentally observed in various physical applications (see [lj|^ and references therein) 
such as radio oscillators, mechanical systems, lasers, electrochemical oscillators, plasma and gas 
discharge, and quantum systems. The study of this phenomenon is also very important from the 
viewpoint of its application to information transmission [^, cryptographic coding with the 
of deterministic chaotic oscillations, and quantum computation [^[^. 

There are several types of synchronization of chaotic oscillations [^: generalized synchro¬ 
nization [^, complete synchronization 10 , antisynchronization 11 , lag synchronization 0 , Re¬ 


use 


quency synchronization |13|, phase synchronization |14|, and time scale synchronization 15 . For 


each type, an appropriate analytic apparatus and diagnostic methods have been developed. Nev¬ 
ertheless, intensive investigations are being continued that are aimed, on the one hand, at the 
examination of different types of synchronization from unique positions and, on the other hand, 
at the search for new types of synchronous behavior that do not fall under the above-mentioned 
types. In spite of the long history of the study of synchronization of chaotic oscillations, many 
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important problems in this field remain unsolved. These include the quantitative analysis of the 
time structure of synchronization of dynamical systems. By this structure we mean spikes in the 
synchronous behavior of the phase variables of systems between which the level of synchronism is 
characterized by a small parameter, i.e., intermittent behavior |16| . 

The concept of intermittency is very important in physics (and not only in physics) for 
the study of structural properties of processes and is not restricted to the synchronism of chaotic 
systems. This fact was noticed even by Mandelbrot (the turbulent flow problem in fluid 
dynamics) and, somewhat later, by Zeldovich and colleagues (problems of chemical kinetics 
and self-excitation of a magnetic field in a random flow of a conducting fluid). For instance, the 
phenomenon of intermittency has application in high-energy particle physics [l^ , cosmology 


19 


m 

the study of rearrangements of attractors in nonlinear dynamical systems 20 - ^ , and in other fields. 
From the physical point of view, intermittency generally implies the emergence of some structures 
of different scales in a medium (for example, vortices, localized strains, temperature irregularities) 
that originally could be absolutely structureless on these scales. From the mathematical viewpoint, 
such behavior is characterized by the presence of rare but strong peaks in the behavior of the 
indicator of a structure - a certain random variable 


16 


The analysis of the structure of synchronism has theoretical significance for nonlinear dy¬ 


namics itself 20 , as well as applied value, for example, in problems of biology and medicine 23 - 25 
stochastic financial mathematics 26 , and so on. The paper |27| on the analysis of the bubbling 


phenomenon is one of the first publications in which the authors experimentally observed intermit¬ 
tency between synchronous and asynchronous behavior in a system of coupled oscillators. However, 
when analyzing this phenomenon, the authors did not investigate the time structure of synchro¬ 
nism as such. Eventually, interest of researchers in the structural phenomena arising during the 
synchronization of chaos grew steadily. At present, the questions of time, space, and spacetime 
structures of synchronism are intensively studied both from theoretical (model) [^-31 and applied 

and power 


standpoints, for example, from the viewpoint of neurophysiology (neurobiology) IMI 
networks 


33 


32 


For all the significance of the problems of quantitative analysis of the time structure of 
synchronization of nonlinear systems, progress in this field falls far behind the investigation of purely 
spatial synchronization patterns in distributed systems. This is clearly illustrated by the above- 
cited papers 23,28-33 . Among the main reasons, in my opinion, are two objective difficulties. 
First, this is the lack of a unified and complete theory of nonlinear dynamical systems pOl. Second, 


during the analysis of spatial patterns, one can extensively apply both the apparatus of network 
theory and the apparatus of the theory of pattern recognition on images (two- and three- 
dimensional scalar fields of physical quantities). In turn, during the analysis of a time or spacetime 
structure, one faces the problem of transition from time series of phase variables to the graphs or 
images that characterize the intermittency of the synchronization process. Not least of all, this is 
associated with the underdevelopment of appropriate tools, because researchers have mainly been 
dealing with the so-called integral coefficient of synchronism, thus studying the synchronization 
parameters on average. Accordingly, the measurement techniques have been developed, as a rule, 
within the same concepts. Note that the paradigm of studying systems and processes on average 
has been widely accustomed in many experimental fields of science. The reasons for the wide 
use of this approach and its fundamental restrictions are presented in detail, for example, in the 
survey 16 . 

In this paper, we develop an original method for the diagnostics and quantitative mea¬ 
surement of the characteristics of intermittency regimes during synchronization of chaotic systems, 
which is aimed at the integrated study of the time structure of synchronism through the analysis of 
the so-called T-synchronization 35 -37 . The method is based on the formalism of symbolic CTQ- 


analysis proposed by the present author [38|[^ (the abbreviation CTQ stands for three alphabets 
with which the method operates: C, T, and Q). One should note that symbolic dynamics, for all 
its seeming external simplicity, is a very strongly substantiated tool for the analysis of nonlinear 
dynamical systems |40| and allows one to investigate complicated phenomena in systems such as 
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chaos, strange attractors, hyperbolicity, structural stability, controllability, etc. (see, for example, 
and references therein). 

In the present study, we propose to consider the relationship between synchronous and 
desynchronous domains integrally. This allows us to carry out a detailed analysis of intermittency 
regimes during synchronization of chaotic systems. Moreover, we essentially revise the principles of 
the symbolic CTQ-analysis: we formulate the encoding rules for the symbols of the base alphabet 
in a rigorous and formal form, which allows us to form a complete set of symbols. In addition, we 
propose a number of measures for evaluating the time structure of synchronization. As an example, 
we take a system of logical mappings, which, on the one hand, is a standard object of nonlinear 
dynamics, and, on the other hand, with regard to the Feigenbaum universality, many results of the 
analysis of this system are extended to a wide class of both model and real objects. 


40-42 


2. Definition of T-synchronization of chaotic systems 


Introduce a trajectory of a dynamical system defined as a discrete sequence (time se¬ 
ries) {sfc}^^, where the phase variable s of the system has dimension N and the trajectory consists 
of K time samples. In this case, each A:th sample can be assigned a time instant tk- 

Define the original mapping that encodes the form of the nth component of the sequence 


in terms of a finite T-alphabet pI8,39 




(n) 




Tt'X, 


( 2 . 1 ) 


Strictly speaking, mapping (2.1) is defined by the relations 
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( 2 . 2 ) 


The graphic diagrams illustrating the geometry of the symbols for the fcth sample 

and nth phase variable are shown in Fig. 

Thus, the T-alphabet includes the set of symbols 


T"^ = {TO, Tl, T2, T3N, T3P, T4N, T4P, T5N, T5P, 


T6S, T6, T6L, T7S, T7, T7L, T8N, T8P}. (2.3) 
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Figure 1. Graphic diagrams illustrating the geometry of the symbols of for the fcth sample and the 

nth phase variable. 


This formula shows that a symbol is encoded as Ti, where i is the right part of 

the symbol codes of the alphabet To*^. In turn, a symbol is encoded in terms of Tii • • • 
see (2.3). The complete alphabet that encodes the form of the trajectory of a multidimen¬ 


sional sequence consists of 17^ symbols in total. 

Now, suppose, for simplicity but without loss of generality, that the time sequence 
of dimension N is formed by the combination of the phase variables of N one-dimensional dynamical 

in) 

systems; i.e., suppose that ' is the value of the phase variable of the nth system at the /sth instant 
of time. 

We will assume that the dynamical systems are synchronous at time instant k in the sense 
of T-synchronization |35| if the condition = 1 is satisfied, where 


* = r = -■ = rru = ■ ■ ■ = 

1 0 otherwise. 


(2.4) 


Taking into account possible antisynchronization |11| between the systems, we should also 
consider all possible variants of inversion of their phase variables: —?• — 1 • In this case, for 

the nth system, a change of symbols of in the /sth sample occurs by the scheme 


TO TO, 

T1 T2, T3N T5P, T3P T5N, T4N T4P, (2.5) 

T6S O T7S, T6 O T7, T6L O T7L, T8N O T8P. 

Denote each variant of inversion by number m. The total number of variants of inversion is M = 

2^“k 

Synchronization between systems can also be set in the lag mode [^. To detect this 
synchronization, one should move a little the phase trajectories of the systems with respect to each 
other ((shifts ^ 0): 

{rrii ^ ..., ir”!™ ^ . (2,6) 


The antisynchronization and lag synchronization modes may coexist; therefore, when cal¬ 
culating a partial integral coefficient of synchronism, we take this fact into consideration: 


5 


S _ 

m,h 


1 

K* + l 


K* 

— Y, Jk\{m, h}, 

k=k* 


(2.7) 


where A:* = 1 -|- max {hi, ..., Hn), K* = K + min {hi, ..., /iat), and K is the length of the se¬ 
quence 

On the basis of the partial coefficient, we calculate the total integral coefficient of synchro¬ 
nism of the systems: 


5^ = max max j^, 

m h ’ 




( 2 . 8 ) 
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i.e., we take a combination of shifts between the trajectories of the systems and a variant of inversion 
of their phase variables that, taken together, provide the maximum number of samples /c, satisfying 
the condition = 1. 


It follows from the definition of the synchronization condition (2.4) that the analyzer pro¬ 
posed evaluates the complete synchronization level 10 and detects antisynchronization [11] with 
lag synchronization (l^ precisely in the alphabetic representation However, according to the 


definition of the geometry of the symbols of the T-alphabet (2.2), a complete synchronization at the 
level of the samples of is a wider phenomenon compared with the complete synchronization at 
the level of - the samples of the sequence itself. The T-synchronism of dynamical systems (with 
respect to the set of phase variables s) is considered from the viewpoint of the shape (geometric 
structure) of the trajectories of the systems in the extended phase space. By the shape (geometric 
structure) of a trajectory of a dynamical system in the extended phase space is meant its certain 
invariant under uniform translations and dilations of the trajectory in the space of phase variables. 
Thus, in a sense, the T-synchronization deals with the topological aspects of synchronization of 
dynamical systems 40 43 . Hence, this opens a potential possibility for the application of the 
analyzer proposed to the study of generalized synchronization of chaos M. 


3. Time structure of synchronization of chaotic systems 


The quantity 5^ introduced in (2.8) characterizes the synchronism of the systems on average 


over a period of tx — ti. As mentioned in the Introduction, most investigations on the synchroniza¬ 
tion of chaos are usually restricted to this situation. However, often a researcher may be interested 
in the time structure of synchronization of systems. Recall that by this structure one means the 
spikes in the synchronous behavior of the phase variables of the systems between which the syn¬ 
chronism level is characterized by a small quantity, i.e., one means intermittent behavior |16| . 

In [^, the present author introduced the concept of a synchronous domain SD - a set of 
samples of a time series that satisfy the condition (V is the symbol of the logical operation OR) 


SDj. : !<//;/ — 1, J^n — 0 V — 0, Jk'" — 0 V — AT -|- ij", 


k' = + Lr, k" = bf^-l, fc'" = + 1 , 


(3.9) 


where 6®, and r are the emergence time, the length, and the ordinal number of a synchronous 
domain, respectively. In this case, the following conditions are satisfied: ^ AT, and the total 

number of synchronous domains (in the original sequence) ^ {K + 1) div2. 

To quantitatively describe the structure of synchronization of systems, the author intro¬ 


duced in 35 a spectral density function of synchronous domains SD: 




[A] = > A], 


SD 


r=l 


(3.10) 


where 5[o, o] is the Kronecker delta and A = 1, A'. 

To analyze the degree of degeneracy of a structure of synchronous domains, we addition¬ 
ally define a quantity - the entropy of the structure of synchronous domains (according to 
Shannon) [^, which makes sense for <5^ > 0: 

^SD = _^pSD i^pSD pSD ^ 

*=1 E M 

i=l 

It follows from Shannons entropy properties that the entropy is minimal (A®^ = 0) when 
the spectrum H^^[L] is degenerate (all synchronous domains have the same length) and maximal 
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in the case of a uniform comb spectrum H^^[L] with a limit number of different 
lengths of synchronization domains equal to 


^cmb = min 


Vl + 86^K-1 


K -6^ K + l 


= inwSD 


crabi 


(3.12) 


where [aj is the integer part of a. 

On the basis of (3.11) and (3.12), we define the relative entropy of the structure of syn¬ 
chronous domains: 

(3.13) 


ASD _ 

i\E — 


^SD 


It makes sense to apply the quantity A|P when the researcher should compare synchronization 
cases that differ in the values of <5® and/or K. 

Nevertheless, for the full description of the intermittent behavior of chaotic systems during 
synchronization, it is obviously insufficient to study only synchronous domains SD. To obtain a 
complete and closed idea of the time structure of synchronism of dynamical systems (a complete 
and closed representation of the intermittency structure), in this paper we additionally introduce 
the concept of a desynchronous domain SD - a set of samples of a time series satisfying the condition 


SD^ : = 0, Jk" = 1 V k" = 0, Jk'" = 1 V k"' = if! -|- ij- , 

k' = ftSD, ftSD + k'' = hf^-l, k'" = hf^ + Lf^ + 1, 


(3.14) 


where 6®'^, and r are the emergence time, the length, and the ordinal number of a desyn¬ 
chronous domain SD , respectively. 

We will schematically illustrate the necessity of a consistent analysis of desynchronous and 
synchronous domains by a model example (see Fig. [^). Let us define three variants (A, B, and C) 
that have the same integral coefficient of synchronism, the same set and sequence of synchronous 
domains SD, but essentially different sets and sequences of desynchronous domains SD. 

If these model variants are compared only with respect to the parameters of SD domains, 
then a wrong conclusion is made on the identity of synchronization processes in these cases. How¬ 
ever, Fig. points to the contrary. Therefore, in addition to SD domains, we consider SD domains. 
Figures Ep and show that, in the case of A, all desynchronous domains have a sample length 
of 2 and, in the case of B, desynchronous domains have a sample length of 2 and 11 (among which 
there is a single desynchronous domain of length 11). The case C has the most complex structure: 
in addition to the fact that desynchronous domains have four different lengths of 1, 2, 3, and 4 
samples, the structure of SD e-)- SD transitions is rather nontrivial. Thus, a consistent analysis 
of the parameters of SD and SD domains, as well as of the structure of SD -H- SD transitions, 
provide much more information on the character of rearrangement of the structure of attractors 
and on the intermittent behavior in synchronized systems. Note that, in most cases (see, for ex- 
), it is synchronization patterns that are subjected to analysis as a rule, while 
the desynchronization areas are considered as a background signal; i.e., they are removed from the 
analysis. 

A quantitative analysis of the time structure of synchronization is possible through the study 
of a graph (see Fig. [^) whose vertices are domains SD and SD with lengths L and L, respectively, 
and arcs (oriented edges) are SD —)• SD and SD —)• SD interdomain transitions. The labeling of 
the vertices of this graph is made by the probabilities and and the labeling of edges, by 
conditional transition probabilities: 


ample, 23 28 - 33 


_E E '] 


L] 

K K I T’ 

E E i] 

1=11=1 


(3.15) 
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Figure 2. A model example: (a) sequence order of desynchronous SD and synchronous SD domains; (b) 
spectral densities of desynchronous and synchronous domains; and (c) scheme of transition SD o SD 
between domains of certain length. The letters A, B, and C denote three cases that have different structure 
of desynchronous behavior. 


where and - are the spectral densities of the interdomain transitions SD —?■ SD and SD —)• 
SD, respectively. 

The analysis of the graph of the SD •(-)> SD transitions opens wide possibilities for the 
application of the currently intensively developing apparatus of network theory 34|44 - 46 , involving 
the methods of statistical physics |47| and symbolic dynamics [40-42 , to the analysis of the time 


structure of synchronization of chaotic oscillations. Note that these topics are the subject of our 
current research. 


4. Analysis of a system of logistic mappings 


A system of coupled logistic mappings is one of popular models of nonlinear dynamics for 

Its theoretical value is associated with the fact that, in 


studying chaotic synchronization [38,48 


spite of its relative simplicity, a logistic mapping gives rise to a wide spectrum of complex, including 
chaotic, oscillation modes [39 49 . A transition to such oscillation modes occurs by the classical 


period-doubling scenario. With regard to the Feigenbaum universality, many results are extended 
to a wide class of both model and real physical, biophysical, chemical, and other systems, which 
also arouses applied interest in the logistic mapping [4^[50| . 

Consider a system of two unidirectionally coupled logistic oscillators: 


Xk+i = 4 Axfc (1 - Xk), 

Vk+i = AX[yk + -f{xk- Vk)] (1 - [yk + 7{xk - Vk)]), 


(4.16) 


where: A G (0, 1] is a control parameter, 7 G [0, 1] is the coupling coefficient of the systems, and 
X, y G (0, 1) are the phase variables of the master and slave systems, respectively. 


The analysis of mappings (4.16) was carried out on the interval fe G [l x 10^, 2 x 10^]. 
Such a shift from k = 1 is explained by the need to reduce the spurious effect of the transient process 
to minimum. Moreover, all the estimates of the quantities were averaged over 1000 realizations 
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of the initial conditions: xi = and yi = ^ 2 , where ^ 1,^2 £ ( 0 , 1 ) are uncorrelated uniformly 
distributed pseudorandom variables. This allowed us to neutralize the memory effect induced by 
the initial conditions on the trajectories of the processes x and y. The value of the coupling 
parameter was varied in the interval 7 G [0, 0.5] with a step of 1 x The value of the control 

parameter was taken to be A = 0.95, which means the existence of a developed chaos regime 
in (4.16) 50 . Note that, according to the results of |39| , the final qualitative complexification 
of the chaotic trajectory occurs at A = 0.98805.... Nevertheless, the choice of the value of the 
control parameter is attributed to the necessity of mutual analysis and coordination of the results 


In particular, it was established in pl8,48 


obtained in the present study with the results of 38,48 
that, for 7 ^ 'jrb — 0.38, a robust regime of complete synchronization is set in the system, and, 
at the point 'ybb — 0.35, a nonrobust regime of complete synchronization is set that persists on the 
interval 'jbb ^ 7 < 7 r 6 - There are four more points on the scale of the coupling parameter that are 
responsible for various rearrangements of the structure of the attractor: 'ywm — 0.0639, 7 ^^ ~ 0.14, 
7 ri ^ 0.2606 and 7^2 — 0.2755. For example, at an unstable quasiperiodic motion is fixed in 
the attractor on which the trajectory stays for a considerable period of time 48 , while, at 7 ^ 1 , the 
trajectory of the slave system (y) contains the maximum concentration of relaxation oscillations, 
which are missing in the master system (x) [38] . 

Figure]^ demonstrates the graph of the integral coefficient of synchronism as a function 
of the coupling parameter 7 . For 7 ^ ~ 0.0191, the lower quantile 6^ becomes greater than the 
upper quantile 6^ for the first time. Thus, in a first approximation, the point 7 ^ can be considered 
as a statistically significant threshold for the emergence of T-synchronization in system (4.16). 
Note that the critical intervals were constructed by the empirical distribution functions of the 
calculated characteristics. The upper and lower quantiles are constructed as two-sided quantiles of 
orders 1 — a/2 and a/2, respectively. The level of statistical significance was taken to be a = 10“^. 
As an estimate for the mean value of the characteristics, we used a median, as a more robust 
indicator compared with the arithmetic mean 51 . Figure demonstrates the graph of A5^ = 
5^ — (5^ as a function of the parameter 7 . 



0.0 0.1 0.2 0.3 0.4 y 0.5 


b 



Figure 3. (a) Integral level of synchronization (the inset shows the region of small values of 7 ) and (b) 
the width of the interval A(5® = 5® — <5®, as a function of the coupling parameter 7 . Red curves represent an 
estimate for the median, and Blue curves, the interval of values with boundaries at probability 1 — a with a 
significance level of a = 10 “^. 


Figure shows that the behavior of the integral coefficient of synchronism (5® allows one 
to reliably detect only a domain of bubbling behavior of the system near the point 'jbb at which a 
nonrobust regime of complete T-synchronization is set. The nonrobustness manifests itself in a more 
than fivefold stepwise broadening of the interval A(5® (see Fig.j^). Much more complete information 
on the rearrangement of the structure of attractors and the time structure of synchronization is 
provided by the analysis of the medians of the spectral densities of the domains SD and SD, which 
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are 


demonstrated in Fig. 
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Figure 4. (a) and (b) as a function of the coupling parameter 7 . The inset shows the function in the 
region of domains of large length ^ 99900. 


The time instant at which complete T-synchronization is set in system ( |4.16 ) at 'jbb can 
also be determined from the analysis of (see Fig. [^). This time instant manifests itself in the 
burst of the lengths of synchronous domains (see the inset in Fig. Ef and Fig. [^), as well as in the 
transition of desynchronous domains into synchronous ones. 
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Figure 5. (a) for domains of length ^ 50 and (b) iJ®^ for domains of length L®^ ^ 80, as a 
function of the coupling parameter 7 . The color intensity scales correspond to those in Figs. |4k and Hb. 


The comparison of Figs and Do shows a significant difference between the behavior 
of desynchronous SD and synchronous SD domains. For instance, the maximum length of syn¬ 
chronous domains on the entire range of variation of 7 is not greater than maxL^^ = 220 samples, 
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whereas maxL®^ = 10^ samples (after setting a robust regime of complete T-synchronization). As 
the coupling increases from 0 to 'jbb, the lengths of synchronous domains monotonically increase, 
and the maximum of a spectral domain is shifted toward longer domains. This indicates that the 
characteristic is insensitive to the rearrangements of the attractor structure, which occur in 
the slave system (y) up to the onset of the nonrobust regime of complete T-synchronization. In 
turn, the analysis of allows one to fix a number of rearrangements of the attractor structure 
in the system (y) that occur under the action of the master system when the coupling varies from 0 
to 7 ft;,. For instance, the regimes at and 7^1 are accompanied by local spikes of the length of 
desynchronous domains, while the regime at 'j^s is characterized by a pronounced comblike struc¬ 
ture of the spectrum of SD domains of small length (see Fig. [^). This kind of spectrum indeed 
corresponds to an unstable quasiperiodic motion where the trajectory of system (4.16) stays for a 


significant period of time 48 


The comparison of Figs. [ 6 }i and also demonstrates a difference in the behavior of the 
degree of degeneracy of the structure of synchronous and desynchronous domains, depending on 7 
- the coupling between the master and slave systems. It follows from the analysis of the figures 
that, almost up to the setting a nonrobust regime of complete T-synchronization, the relative 
entropy of the structure of domains SD varies very little, remaining near the level of 0.5. There 
are local statistically significant spikes near the rearrangement areas of the attractor structure 'jwm 
and 7 ri. The width of the interval of values with boundaries at probability 1—a with the significance 
level a = 10“^ also remains virtually constant. As the parameter 7 approaches the critical value •jbb 
and the system passes to a nonrobust synchronization regime, the interval of values of A|P with 
boundaries at probability of 1 — a opens up, and the median A|P experiences a jump from 0.5 
to 0.73. A further increase in the coupling force to 7 = 0.5, the width of the interval is not 
changed, while the median A|P decreases to values close to zero. In turn, as the coupling parameter 



Figure 6. (a) A|P and (b) A|P as a function of the coupling parameter 7 . Red curves represent an estimate 
for the median, and Blue curves, the interval of values with boundaries at probability 1 — a with a significance 
level of a = 10“^. 


increases from 7 = 0, the relative entropy of the structure of synchronous domains SD first locally 
(statistically significantly) increases from 0.35 to 0.375 near the value of 'y^m and then returns to 
the original level, and, starting from about 7 ^^, again starts to increase. Before the 7 riregion, there 
is a local maximum, and, in the region [ 7 ^ 1 , 7 r 2 ]) there is a local minimum. As the coupling force 
increases further, the interval of values of A|P with boundaries at probability of 1 — a opens up 
near 7 ftft, and the median A|P, upon reaching a maximum of 0.57, decreases to values close to zero. 

The analysis carried out implies that the relative entropies A|P and A|P are a very sen¬ 
sitive indicator of rearrangements of the attractor structure that occur in the regime of chaos 
synchronization. 
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5. Conclusions 


We have proposed a new method for integrated analysis of the time structure of synchroniza¬ 
tion in nonlinear chaotic systems. The method allows one to diagnose and quantitatively evaluate 
the intermittency characteristics during synchronization of chaotic oscillations in the so-called T- 
synchronization mode |35l|36|. The approach is based on the formalism of symbolic CTQ-analysis 


proposed by the author in 138, 39 


To demonstrate the main capabilities of the method, we have carried out an analysis of a 
system of two identical logistic oscillators with unidirectional coupling that are in the developed 
chaos regime. This system is a standard model of nonlinear dynamics. In addition, in view of 
the Feigenbaum universality, many results of the analysis can be extended to a wide class of both 
theoretical and real objects (49,50 . The analysis carried out has revealed the setting of nonrobust 


and robust regimes of complete T-synchronization in the system. The nonrobust regime is accom¬ 
panied by a bubbling attractor that exhibits intermittency between synchronous and asynchronous 
behavior of the oscillators. In addition, we have confirmed the results of 48 concerning the fact 


that, for certain values of the coupling force between subsystems, an unstable quasiperiodic mo¬ 
tion is detected in the attractor on which the trajectory stays for a considerable period of time. 
We have also established that the regime, revealed in 38 , in which the trajectory of the slave 


subsystem contains the maximum concentration of relaxation oscillations, which are missing in the 
master subsystem, leads to a local increase in the deviation of the length of desynchronous domains. 
This fact increases the intermittency between the synchronous and asynchronous behavior of the 
oscillators. 


The analysis has confirmed that the widely used method in which only synchronization 
patterns are analyzed, while the desynchronization areas are considered as a background signal 
and are actually not analyzed, should be regarded as methodologically incomplete due to the loss 
of important information on the character of rearrangement of the structure of attractors and the 
intermittent behavior in synchronizing systems. 


In conclusion, note that the method considered, which is based on the analysis of T- 
synchronization, can be successfully applied to the study of multidimensional systems consisting 
of two or a greater number of coupled nonidentical oscillators, including multidimensional lattices 
of oscillators with arbitrary topology. The approach described can be applied to the analysis of 
experimental data, because it does not require any a priori knowledge of a system under study. 
Moreover, the invariance of the analyzer with respect to the shifts and dilations of phase trajec¬ 
tories (^ allows one to investigate the synchronization of strongly nonstationary systems. It is 
quite possible that this fact will allow one to effectively apply the method proposed to the analysis 
of multidimensional time series generated by physical and engineering 
financial |26l|37l|46 


biophysical 23 - 25 


power 


33 , and other systems. 
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